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Conversion of entangled states under (Stochastic) Local Operations and Classical Communication 
admits the phenomenon of catalysis. Here we explore the possibility of a copy of the initial state 
itself to perform as a catalyst, which we call a self-catalytic process. We show explicit examples of 
self-catalysis. Necessary and sufficient conditions for the phenomenon to take place are discussed. 
We numerically estimate how frequent it is and we show that increasing the number of copies used 
as catalyst can increases the probability of conversion, but do not make the process deterministic. 
By the end we conjecture that under LOCC the probability of finding a self-catalytic reaction do 
not increase monotonically with the dimensions whereas SLOCC, does increase. 
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I. INTRODUCTION 


The conversion between bipartite or multipartite quantum states through local operations is a central concept 
of entanglement theory. For instance, it is the criterion used to classify entanglement, namely, two states have 
equivalent entanglement if they can be converted to each other [T]. Moreover, it is usual and natural to consider 
a state more entangled than other when the first can access the former by the allowed transformations [2]. Two 
important sets of allowed transformations are the deterministic local operations and classical communication (LOCC) 
and its stochastic version (SLOCC), where the transformation only needs to succeed with positive probability [3]. 
Although this hierarchisation is relatively simple for bipartite pure states 0 iil> it is quite involved if mixed or 
multipartite states are considered ISHZl. 

The problem of convertibility of bipartite pure states has been solved by Nielsen[2], using the concept of majorisation 
|8]. However, Jonathan and Plenio discovered a surprising effect |9]. They have shown the existence of pairs of states 
which are not directly inter-convertible but such that their conversion is possible if another (necessarily entangled) 
state is attached to them. That is, the pair of states have incomparable entanglement but they may become ordered 
if an extra system is attached to the original ones. Such extra state that makes a transformation possible, without 
being consumed, is called a catalyst. More recently, the problem of convertibility has received experimental attention 
m and has also been connected to basic results in thermodynamics El and phase transitions m 

In this paper we explore the following question: can one of the states of a non inter-convertible pair be used as 
a catalyst? We answer this question positively, providing explicit examples. We also explore how common such 
processes are for low dimensional systems. An interesting situation is when a state |T) is not able of self-catalysing 
a transformation, but a number of copies, is. We exhibit examples where more than one copy is required and 

study how augmenting the number of copies can increase the probability of conversion. 

In section |TI| we review basic notions of state conversion and catalysis under LOCC. In section |TI^ we show explicit 
examples of self-catalytic processes, explore, through numerical analysis, how frequent they are under random samples 
of incomparable pairs and how it depends on the size of the systems. Section [TV] review the notions of probabilistic 
catalysis under SLOCC, while the natural questions of self-catalysis under SLOCC are discussed in section [V[ We 
close the main text with final remarks and further problems in section im 
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II. CATALISYS 

We say that a bipartite state \a) G Ha ^ Hb access a state |/3) G (g) if there is some LOCC operation, rep¬ 
resented by a completely positive trace preserving (CPTP) map A, such that A{\a) (a|) = |/3) (;d|, where Ha^H'a^Hb 
and H'q are finite dimensional Hilbert spaces and A maps density operators acting on Ha Hb to density operators 
acting on H'j^ (g) H'^. In such case, we write \a) \/3). If there is no LOCC operation able to convert \a) to |;d), we 

shall write \a) ^ \^). 

For example, for a pair of qubits, the Bell state |$+) = (|00) + |ll))/\/2 can access any two-qubit pure state. Indeed, 
naming our qubits A and H, and writing the target state as \/3) = a |00) -f-6 |11), we can, for example, make A unitarily 
interact with an auxiliary qubit A', so that lOy^O^') a lOy^O^')+ ^ |1 a1a') and i-G^ a + ^ |1 a0a')- Then, 

[(IO^Ob) + \1a1b))/V2] |0^.) ^ (l/\/2)(a |0^0b) + 6 |U1b)) |0^.) + (l/\/2)(a |0^1b) + b |U0b)) |U'). The lab with 
qubit A can make a measurement on the computational basis of the auxiliary A' and send the result to the lab 
holding qubit B. If the result is 0, they already share the desired state, while the result being 1, a NOT operation, 
|0 b) |0b), can be applied to qubit B to also leave the system AB in the desired state. This result 

d-l 

|ii), which justifies |^j") to be called a maximally entangled 

i=0 

state. 

Nielsen, in the seminal paper [2], provided a simple necessary and sufficient condition for determining whether a 
general bipartite state \a) can access a state \/3) in terms of their corresponding Schmidt vectors [T]: 

Theorem 1 (Nielsen Criterion). If a = (oi,...,On) and [3 = (/3i,..., are the (non-increasing) ordered Schmidt 
vectors of \a) and \P), respectively, we have jo) ^ \P) if, and only if, 

( 1 ) 

1=1 1=1 

for 1 < k < min{n,m}. 

By an ordered Schmidt vector A = (Ai,..., A^), we mean that Ai > A 2 > ... > A^. In this work, we shall always 
assume that Schmidt vectors are ordered. When vectors a and [3 satisfy conditions 0 we say that a is majorised by 
(3 and write a < 

Revisiting the example above, we have (1/2,1/2) for the ordered Schmidt vector of and (|ap, |6p) for a |00) -f 
6 111), assuming \a\ > \b\, it is straightforward to apply the criterion and verify that |^+) ^ [a |00) + 6 |11)]. 

A consequence of the criterion is the existence of pair of states \a) and \^) such that \a) \^) and \/3) \a). 

Actually, the only case where such an order is total is for two qubits. For instance, consider a state of two four-level 
systems with Schmidt vector (0.4, 0.4, 0.1, 0.1) and a two-qutrit state with Schmidt vector (0.5,0.25,0.25). Indeed, 

0.4 = ai<Pi= 0.5, (2a) 

0.4 “I" 0.4 = OL\ “h 0 L 2 ^ (3i T {32 — 0.5 -1- 0.25. (^t)) 

In Ref. [9], the authors surprisingly show that it is possible to circumvent such non-accessibility between these states 
by making the parts to share an entangled state |/^) such that ja) (8) |/^) ^ \P) ^\k,) (see Figure Q. Since the state |/^) 
allows for a previously forbidden conversion, but at the end of the process it remains unaltered, it is called a catalyst. 
In this sense, we say that \a) n-access \p) when \a) ^ \/3), but |a) 0 |/^) ^ \^) 0 \n). In this specific example, the 
two-qubit state \k) with Schmidt vector ^ = (0.6, 0.4) is a catalyst. 

Nielsen’s criterion ensures that all information about the (possibility of) conversion is contained in the Schmidt 
vectors. Therefore, we can explore our knowledge regarding probability vectors to provide more examples of catalysts. 
lfa = (0.5,0.4,0.05,0.05), P = (0.7,0.15,0.15), kS = (0.7,0.3), and ^2 = (0.75,0.25), we obtain: 

a ^ P (3a) 

ka ^ /3 (8) /7i (3b) 

d 0 k2 ^ 13 k2. (3c) 

This is a good example of non-unicity of catalysts. Indeed, for a given forbidden transition, \a) \/3), and fixed local 

dimensions for the catalysts, the set of allowed vectors ^ is a polytope m- 

Another interesting geometric fact is that catalysis is possible for every bipartite scenario, starting from 4x3, z.e., 
for all effective dimensions m > n, m > 4, and n > 3 it is possible to choose d = {ai, ^ {Pi,...,/3n) = P 


generalises for a pair of qudits in the state 1 
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FIG. 1: Although the conversion \a) \/3) is not allowed under LOGO, the state |/^) can be used (but not 

consumed) to make it viable. 


a 

/3 

ttCopies (A) 

(0.900,0.081,0.010,0.009) 

(0.950,0.030,0.020) 

N = 1 

(0.900,0.088,0.006,0.006) 


A = 2 

(0.908,0.080,0.006,0.006) 


A = 3 

(0.918,0.070,0.006,0.006) 

(0.950,0.030,0.020) 

A = 4 

(0.925,0.063,0.006,0.006) 


A = 5 

(0.928,0.060,0.006,0.006) 


A = 6 

(0.900,0.081,0.010,0.009) 

(0.950,0.030,0.019,0.001) 

A = 1 


TABLE I: Number N of copies required to make a deterministic catalyst for the process \a) \/3). 


with catalyst k. The essential step, after the previously described examples, is that given a, yd, and we generically 
can increase dimensions by one and construct a forbidden transition a' ^ (3' with the same catalyst by using 
a' = (q^i, ...,(a^ — e, e) and P' = (ydi, ...,yd^ — e',e') with small enough e > 0 and e' > 0. This allow us to construct 
examples for all such m and n. 


III. SELF-CATALYSIS 

In this section we address the main question of this paper for the case of LOCC convertibility: can a bipartite 
quantum state be itself the catalyst of a forbidden conversion? To be more precise, is there a forbidden conversion 
I (a) ^ I yd), such that |(a)G)|(a) ^ \/3)0\a)? Refining a little bit more, is it possible that we still have | a) G)| a) ^ \/3)0\a), 
but a larger number of copies of \a) would do the job, z.e., \a) G ^ |yd) G for some N > 17 

In Table |T] we list examples to answer affirmatively these questions. The first is an example of self-catalysis from 
a two-ququart state to a two-qutrit state. Then, for the same scenario, there is a list of multi-copy self-catalysis, 
with the corresponding minimal number of copies. At the last row, we come back to single-copy self-catalysis from a 
two-ququart state to another two-ququart state. 

Remark. When ja) has 4 non-zero Schmidt coefficients and |yd) has 3 non-zero Schmidt coefficients, the phenomenon 
of self-catalysis can be noted, that is, even in the minimal dimensions to occur catalysis (see m) the self-catalysis 
can also carry out. Geometrically, this mean that even is the smallest dimension (most restrictive) case, there are 
cases where the source state a belongs to the polytope of catalysts for the reaction a ^ [3. Moreover, it is possible 
to construct examples for any scenario through a similar argument presented at the Section |TI| 


A. Stability under small perturbations 

It is important to mention that the phenomenon of self-catalysis, as it happens with catalysis, is generically robust 
against small perturbations of the state vectors involved. That is, suppose that one is aiming to perform a self-catalytic 
process ja) G |(a) ^ |yd) G ja), but it actually implements states ja'), |yd'), where ja') ~ ja) and |yd') ~ |yd). Our claim 
is: generically, if \a) a-access \j3), then ja') a'-access yd'. It is easy to see that this will be true, depending only on 
the inequalities implying that |(a) ^ |yd), as well as those assuring that ja) G |(a) ^ |yd) G |(a), all be strict (except 
the last one, which is granted by normalisation). Denote by |A| some norm {e.g., Euclidean) of the vector A and n,m 
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the sizes of Schmidt vectors a,/?. Now, assuming that we have, for some /, for 1 < /c < / and 

^ Sj=i for / < /c < n — 1, the same set of inequalities will hold for the entries of vectors a' and P' if 

\d — a'\^\P — P'\ < e, as long as e is small enough. This implies that ja') and \P') are still incomparable. A similar 
reasoning can be applied when the incomparability of vectors a and (3 is due to two or more changes of signs in the 
inequalities. In the same way, ^ < nm —1, imply 

if \d — d'\^ \P — P'\ < e, for small enough e, which proves the claim. Naturally, the argument includes the well-motivated 
situation when \P') = |/d) as a special case. 


B. Self-Catalysis under LOCC for random Schmidt vectors 

We have numerically investigated how usual the phenomenon of self-catalysis among pairs of incomparable bipartite 
states is. Fixing the sizes of a and yd, we randomly sample pairs of such vectors until we find incomparable ones. The 
sampling of each vector is done by uniformly sorting unitary vectors in (8) i.e. sorting according to the Haar 
Measure in the respective state spaces Emu], and then calculating the correspondent Schmidt vector. After finding 
a pair of incomparable states, we test whether the first of the vectors can be used as a catalyst for the conversion. 
The results show that for this method of sampling and for small dimensional systems, the phenomenon is actually 
atypical Moreover, the numerical estimations seem to imply that the phenomenon of self-catalysis is atypical in any 
dimension. 



FIG. 2: Probability of finding a pair of states exhibiting self-catalysis as function of the dimension of each system, 
that is P(|(a) G |(a) ^ I/d) G |(a) | \a) ^ \P)). Each symbol means an average over a distinct set of random choices. 

For each size explored there are three different symbols, which indicates a reasonable stability in this sampling 

process. 

Figure[^shows a numerical estimation for ¥{\a) (g) \a) |yd) (g) \a) \ \a) ^ |/d)), that is, the conditional probability 
of self-catalysis given a pair of incomparable Schmidt vectors, as a function of their sizes. Note that this probability 
increases with the size of the Schmidt vectors, until sizes about 20 and then starts to slowly decrease, but in fact for 
any dimension its order of magnitude shows that the phenomenon is present, but rare. 

We believe the global maximum present in the Figurecan be explained by the algebraic character of the comparison 
between vectors and the concentration of measure phenomenon. In order to ja) (g) ja) access ja) (g) |yd), all the 
corresponding inequalities of Eq. 0 must be satisfied. Then, in one hand, by increasing the dimension of the vectors, 
this comparison become harder to be valid (more inequalities have to be obeyed). On the other hand, it is known 
that a measure concentration phenomenon takes place for increasing dimensionality, in the sense that the Schmidt 
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FIG. 3: Mean normalized entropy for Schmidt vectors sorted by Haar measure. Here we have explored the behavior 

of the entropy until effective dimension 100. 


vectors become typically closer and closer to the constant vector This phenomenon can be highlighted by 

the average normalized entropy of the sorted vectors: if this average is close to 1, it means that most vectors are close 

to the constant vector. This average value is know to be exactly i^[Ylk=d-\-i \ ~ conjectured by Page [19] 

and latter proved by Foong and Kanno m 

If the vectors of the pair are incomparable, that is, some inequalities of Eq. Q are not satisfied, but both are 
close to the constant vector, it will be easier for the catalyst to make the transition possible, since the corresponding 
inequalities for the vectors with the catalyst attached will be easier to be satisfied. 

We can note that the concentration of measure increases very fast at low dimensions (see Figure]^, helping the 
possibility of “organising” the Schmidt vectors with the catalyst attached, which justifies the corresponding increasing 
in the probability for self-catalysis. However, around dimension 20 the concentration happens much slower and 
presumably is not fast enough to overcome the size effect (which makes majorisation more difficult), so the probability 
of finding a self-catalytic pair decreases after this point. 

Summing up, the probability seems to be a smooth function of the dimension, reaching its higher value for dimensions 
of each Schmidt vector close to 20 and, apparently, converging to a value considerably smaller than 1. 

Remark 1. We have restricted the analysis above to events where ja) ^ \P). For large dimensions we have checked 
that this event is typical. Indeed, in such regime the entries of vectors a and before ordering, are essentially 
concentrated random variables fluctuating around a fixed value. If we look, then, to Eq. 0 we see that there is a 
good chance for the sign of the inequality to change as we vary the index k, implying that the states are incomparable. 
Therefore, we expect that P(|Qf) ^ |;d)) ^ 1, for 1. 


IV. CATALYSIS UNDER SLOCC 

It is possible to generalize the concept of accessibility if we allow for non-deterministic processes. In this case, we 
can look for the probability Pslocc{\<^) \P))^ or P^d^) ^ \fi)) for short, of having the state conversion ja) ^ \fi) 
under the best local strategy, i.e., optimizing P under the conditions defining SLOCC. It is interesting to recall 
that the famous result on inconvertibility between W and GHZ states refers to such conditions on the multipartite 
scenario m 

As shown by Vidal m, for the bipartite case, the Schmidt vectors also encode this maximal probability of conver¬ 
sion, Ps{\(^) 1/^)) through 

Theorem 2. Let a = (oi,..., and be ordered Schmidt vectors for states jo) and \/3), assuming 
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k-l 

> 0- Define Ek{X) = 1 — Then, the optimal transformation probability is given by 

1=1 

= (4) 

For instance, if 

a = (0.6, 0.2, 0.2) and fi = (0.5,0.4,0.1) (5) 

are the Schmidt vectors for the two-qutrit states ja) and \fi), respectively, we get: 

Ps{|a) ^ 1/3)) = 0.8, (6a) 

Ps(|/3) ^ |a)) = 0.5. (6b) 

The following Proposition shows that the optimal probability of conversion Ps attains 1 precisely when ja) access \fi). 

Proposition 1. Let a and /3 be a pair of random independent Schmidt vectors with same size n, then the event 
{|o) ^ \/3)} is equal to event {P^do) ^ |/3)) = 1}. In particular P(|(a) ^ |/3)) = P(P 5 '(|(a) ^ |/3)) = 1). 

Proof Suppose that \a) \/3), thus: 


'ik e {1,2, ...,n}. 


i=l 


i=l 


Then Ek{a) = 1 
therefore Ps{\(^) - 


Et'Pfe > 1 


EtPfe = EkW),yk 


\f3)) = 1. The converse is similar. 


G {l,2,...,n — 1} with the equality holding if k 


( 7 ) 


1 , 

□ 


Note that from our considerations at the end of section IIIB we expect that P(|(a) ^ \f3)) ~ 0 (using LOCC) for 
large n, since the event {|o) ^ \/3)} is in the complement of {\a) ^ \/3)}- 

We numerically estimate the SLOCC average rate of conversion between random states. For this, we generate 
incomparable Schmidt vectors following the Haar measure and calculate the probabilities of conversion from the first 
to the second and also the maximum conversion probability. The results shown in Figure [^clearly shows that, given 
a random pair a, it is very common to have a large probability of conversion from some of them to the other, that 
is E[max{P 5 (|a) ^ \fi) ,Ps{\fi) |<^))}] ^ 0.8. Moreover we have a smaller, but still significant, average probability 
of conversion E[P 5 '(|Qf) ^ slightly below 0.6. 

Also in the probabilistic scenario the presence of an extra state can improve the probability of conversion between 
two states m- As a chemical catalyst, this extra state is used, but not consumed, to increase the rate (probability) 
of a reaction (conversion). In the above example, if k = (0.65,0.35), we arrive at; 


a(g>K = (0.39,0.21,0.13,0.13,0.07,0.07), 

/3 0 K = (0.325,0.26,0.175,0.14,0.065,0.035). 


(8a) 

(8b) 


Therefore P^da) (g) |/^) ^ |/3) (g) |/^)) 0.904, and |/^) can be viewed as a probabilistic-catalyst in the stochastic scenario 

for the conversion that starts in a and ends in despite the fact that Ps{\fi) ^ \f^) \<^) ^ |/^)) = 0-5, and then |/i:) 

does not increases the probability of conversion for the transformation that starts in and ends in a. In this sense 
Jonathan and Plenio pointed out [9] that if Ps{\(^) under the best local strategy, is equal to an/fin^ then this 

probability can not be increased by the presence of any catalyst state. Feng et. al m improved this result, obtaining 
the following theorem: 

Theorem 3. Let a and be two n-dimensional probability vectors written in non-increasing order. There is a 
probability vector R such that Ps{ct 0 H ^ 0 il) > Ps{ct ^ if, and only if. 


Ps{a -)> /3) < min { -E, 1 

L Pn 


( 9 ) 


A deeper connection between LOCC catalysis and its stochastic counterpart connects the probability of occurrence 
of the event {o ^ fi} and the maximal probability of stochastic conversion, Ps{\c^) C |0) ^ \I3) C 




7 



FIG. 4: Sampled probability of conversion under SLOCC for randomly chosen (following the Haar measure) 
incomparable states as function of the dimension of each system. Each symbol represents an average over a distinct 
set of randomly chosen pairs. There are two sets of three symbols for each explored dimension. For each dimension, 
the smaller results consider conversion from the first to the second E[P 5 '(|a) ^ |/^))], and the larger the maximum 

conversion rate E[max{P 5 (|(a) ^ \P) ^Ps{\P) |<a))}]. 


Proposition 2 . Let a and ^ be a pair of random independent Sehmidt veetors with same size n. Then 


Proof. 


“[suplP^da) O 10 ) ^ 1 / 3 ) O 10 ))} > Ps'(la) |/?))] > ^ - IP(|a) I/?))- 


P[P5(|a) 0 10) ^ 1/3) 0 10)) > Ps{\a) ^ 1^))] = 1 - P{Ps{|a) ^ |^)) = minK//3„, 1)} 

> 1 - P{Ps{|a) ^ 1^)) = an/M - P{|«) ^ 1^)1 

1 


>__P{|a)^|^)}. 


( 10 ) 


( ll a) 

( ll b) 

( ll c) 


Where (fTTal 


( |lla| ) comes from Thm i H from set theory and Propositionand finally ( |llc| ) from {P^do) ^ fd)) = 
C^n/M ^ l^n < fdn}- □ 

Remark 2 . From Remark^ we know that as n grows, E(|(a) ^ \P)) ~ 0, so we eonelude from Proposition^ that 
P[Ps((S^/3)<min{|^,l}]>l/2. 

In Ref. [17] a necessary and sufficient condition for a state \k) to work as a probabilistic catalyst was provided: 

Theorem 4. Suppose that a and fd are two non-inereasingly ordered n-dimensional probability veetors, and P{d 
/3) < min 1^, l|. Define 


L = I l;l < I < n, and Pia -^ ^) = I . ( 12 ) 

I 

Then a non-inereasingly ordered k-dimensional probability veetor k serves as a probabilistie eatalyst for the eonversion 
from jo) to \/d) if, and only if, for all ri, r 2 ,..., G P U {n + 1} satisfying ri > r 2 > ... > 7 ^ n + 1, there exist i and 







j, with I < j < i < k, such that 


Ki Prj 

- < - — 

— 1 


or 


> 


Pr,-1 


(13) 


By definition, whenever one of the inequalities (T^ includes an index n 1, it is considered to he violated, so the 
other one must necessarily he satisfied. 


We should stress a couple of facts about the set L: first of all, note that it is a key ingredient for identifying catalysts 
for a given conversion, since in a certain sense it determines which indexes are really important for the comparison 
between and n. Secondly observe that typically L has only one element /, i.e., the minimum which determines the 
probability of conversion (see Thm. 2) is non-degenerate. 


V. SELF-CATALYSIS UNDER SLOCC 

The phenomenon of self-catalysis can also take place when considering conversions under SLOCC. Namely, if a 
conversion \a) \/3) takes place with optimal probability 0 < p < 1, it can be the case that the optimal probability 

for I (a) (8) I (a) ^ |/3) O |(a) be p' > p. Indeed, for the same Schmidt vectors a and given by Eq. (§, there is a gain in 
the probability of conversion if we use the state \a) itself as a catalyst. Using Eq. Q, we have: 


Ps{\a) (g) |a) ^ \fi) (g) \a)) - 0.889 > 0.800 = Ps{\a) 


(14) 


As it happens in the context of LOCC operations, the conversion between states can depend on the number of 
attached copies of \a). Table [H] shows, for the example we are considering (given by Eq. (§), how the probability of 
conversion increases with the number of copies of |(a) to be used as catalyst. 

This example may suggest that, by increasing the number of copies of ja), the probability of conversion approaches 
one. This is not always the case, however. Note that we can bound from above the probability in Thm. for the 
pair \a) (g) and \I3) (g) since Ps{\(^) ^ ^ \fi) 0 < an/fin^ for all N > 1. Therefore, as long 

diS an/fin < 1, no matter how many copies of \a) we have, the probability of conversion will not exceed an/fin- The 
previous reasoning allows us to state the following Proposition: 

Proposition 3. Let a and fi he a pair of ordered Schmidt vectors with n non-null components. If an/ fin < 1 ; then 


Ps{\(y) (g) |o) 




^)®iV) < V A^>o. 

fin 


a 


(15) 


Eor example, given a = (0.60,0.21,0.10,0.09) and fi = (0.55,0.25,0.10,0.10), we have a pair of incomparable 
states with Ps{\(^) |/5)) = 0.88, but since a^/fi^ = 0.9, the probability of conversion under SLOCC using self¬ 

catalysis is limited by 0.9 and indeed, for this case, A = 1 is already optimal, since PsiM C) \ a) \fi) C |o)) = 0.9. 
Analogously, for the pair a' = (0.40,0.34,0.15,0.11) and fi' = (0.50,0.21,0.17,0.12), Table III shows the behavior of 
Ps{W)^W)^^ \ fi')^W)'^^) with the number N of copies of \a'), and since a!^/fi'^ =_^1667 < 1 the probability 


of conversion may increase, but can not reach 1. Moreover, there is the case (see Table IV) where, for a given pair 
a, fi of incomparable Schmidt vectors, the probability of conversion increases monotonically with respect to N and, 
for an Aq < oo, Ps{\(^) C) ^ \fi) (g) = 1 and, by Proposition \a) (g) ^ \fi) (g) 

A particular case of Theorem allows us to obtain a necessary and sufficient condition for having probabilistic 
self-catalysis for a single copy: 

Criterion Let a and fi he two n-dimensional Schmidt vectors with Ps{ct fi) < min l| and 


/; 1 < / < n, and Ps{^ fi) 


Ei{a) 

Efifi) 


(16) 


The vector a serves as a prohahilistic self-catalyst for the transformation from jo) to \fi) if, and only if, for all 
ri, r 2 ,..., G 1/ U {n -h 1} satisfying ri > r 2 > 7 ^ n + 1, there exist i and j, with 1 < j < i < n, such that 


fir 


^ < 

fi^i — 1 
or 


(17) 


(18) 


firj-1 

- > — - - . 

aj fir- 
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tt Copies (N) 

Ps(|a) ® |a)®^ ^ 1/3) ® |a)®^) 

0 

- 0.800 

1 

- 0.889 

2 

- 0.907 

3 

- 0.926 

4 

- 0.932 

5 

- 0.940 

6 

- 0.944 

7 

- 0.947 

8 

- 0.950 

9 

- 0.952 

10 

0.955 


TABLE II: Increase of optimal probability for the conversion \a) \/3), where a = (0.6,0.2,0.2) and 

(3 = (0.5, 0.4, 0.1), with the number of copies N of ja) used as a catalyst. 


tt Copies (N) 

Ps{\a') ® ^ 1^') l®')®"^) 

0 

0.8965 

1 

0.9038 

2 

0.9072 

3 

0.9092 

4 

0.9105 

5 

0.9109 

6 

- 0.9110 


TABLE III: Increase of optimal probability for the conversion ja') ^ \P')^ where a' = (0.40,0.34,0.15,0.11) and 
P = (0.50,0.21,0.17,0.12), with the number of copies N of ja') used as a catalyst. 


Whenever one of the inequalities 0 or has an index n + 1, it is eonsidered to he violated, so the other one must 
necessarily he satisfied. 


Remark 3. Again, we observe that typically L has only one element I, since the minimum which determines the 
prohahility of conversion is non-degenerate with prohahility 1. 


tt Copies (N) 

Ps(|a)®|a)®^ ^ 1/3)® |a)®^) 

0 

- 0.600 

1 

- 0.818 

2 

- 0.911 

3 

- 0.957 

4 

- 0.981 

5 

- 0.994 

6 

= 1 


TABLE IV: Increase of optimal probability for the conversion \a) \P), where a = (0.928,0.060,0.006,0.006) and 

P = (0.950,0.030,0.0195,0.0005), with the number of copies N of \a) used as a catalyst. 
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A. Self-Catalysis under SLOCC for random Schmidt vectors 


The criterion above, together with Proposition and the behaviour of L, have interesting consequences for the 
probability of occurrence of self-catalysis. 

From Remark 1^ we know that the event [P^d^) ^ \P)) < min{^,l}] has probability > 1/2. Conditioning on 
this event we can then analyze the validity of Ineqs. ( pT| ) and ( p!^ . With probability 1 we must have L = {/}, for 
some 1 < / < n. Therefore, we can choose the indexes in only two ways: either ri = r 2 = ... = Vn = / or 


n + 1 and = ... 


/, for some j. For the second case, one can always satisfy one of the 


ri = r2 = ... = ' 

inequalities using the index n +1. For the first case, the r.h.s of the inequalities always have index I and we can lower 
bound the probability for at least one of Inequalities (18) to be valid: 


Qfi /3/_i 
max I — I > 


c^j 


Pi 


> P[ max 

l<j<n OLj 


■-}> rn^{^}] = 1/2, 

l<j<n Pj 


(19) 


using that maxi<j<^<^{^} = maxi<j<^{^^^}, since a is ordered, and the fact that the two random variables on 
the second term are independent and identically distributed. Putting these together we get that the probability for 
having SLOCC self-catalysis is > 



FIG. 5: Probability of finding a pair of states exhibiting self-catalysis under SLOCC, 

P[P 5 '(|(a) 0\a) \l3) 0\a)) > P(|a) ^ |/^))], as function of the dimension of each system. Each symbol represents an 

average over a new set of randomly (Haar) chosen pairs. There are three symbols for each explored dimension. 

We have also numerically investigated the typicality of probabilistic self-catalysis by 1) sorting a pair of incomparable 
Schmidt vectors of same fixed dimension; 2) counting how many of them show the effect; and 3) computing the average 
gain in probability. To be more specific, we consider as a success case the situation where the pairs are such that 
Pi = P^d^) ^ I/d)) < P^d^) C |<a) ^ I/d) C |<a)) = P2 and we compute the average value of p2 — Pi- In order 
to avoid counting cases where pi < p 2 due to numerical fluctuations, we only consider as valid those vectors where 
P 2 > (l + 10“^) pi. The results are shown in Figures]^ and 

First note that they are consisted with the lower bound oil/4 estimated before. Comparing with the deterministic 
case, probabilistic self-catalysis is much more frequent, as expected. Even more, here we do not have the same qualita¬ 
tive behaviour: the probability of having a pair of incomparable states exhibiting self-catalysis increases monotonically 
with the size of the Schmidt vectors and seems to be converging to a value around 0.6. Meanwhile, this is not the 
behaviour of the average probability gain, shown in Eigurej^ The average gain in probability is relatively small for all 
system sizes, and decreases even more for larger sizes. But note that we consider only one copy of \a) attached, and 
while for many copies the average gain must be larger. It is important to recall Eig. however, which tell us that a 
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FIG. 6: Average gain of probability, E[P 5 '(|(a) G \a) \^) G |(a)) — Ps{\c^) IP))]^ as a function of the size of 

Schmidt vectors (sampled following Haar measure), considered only those pairs where self-catalysis occurs, z.e., 

Ps{\(^) G |tr) ^ 1^) G \a)) > Ps{\c^) \P))- 


pair Pj has, on average, a probability of direct conversion close to 0.6, which naturally bounds the catalytic gain 
to about 0.4. 

Finally, Figure represents, for randomly chosen pairs of incomparable Schmidt vectors with size n = 45, the 
self-catalytic probability gain versus direct conversion rates. The diagonal straight line just represents saturation of 
probability. Some concentration close to the horizontal axis is natural, representing the cases where self-catalysis does 
not happen. However it is not clear why there is the bold concentrated cloud in red, where the majority of the pairs 
fit. It empirically means that the most typical situation for a pair of Schmidt vectors of the same size is to have a 
large probability of conversion and to have a considerable (but not maximal) self-catalytical gain. 


VI. CONCLUSIONS 

In this paper we have shown the possibility of self-catalytic entanglement conversion for both LOCC and SLOCC 
scenarios by providing explicit examples of them. We have explored numerically and by arguments of typicality 
how frequent they are by showing that it is much more common in SLOCC case than the deterministic one, but 
despite the fact that self-catalysis under SLOCC become more common as the systems sizes increase, the direct 
self-catalysis decreases with the dimension. Moreover, we also investigated how the phenomenon may depend on the 
number of copies used as catalyst, finding examples of different behaviours, running from cases where there is no 
gain in considering multiple copies, to cases where the conversion becomes deterministic for a finite number of copies. 
Since we rest on numerical techniques, we could not guarantee the existence of a transition where the probability 
of conversion asymptotically goes to I. Finally, we computed the average gain in probability in the SLOCC case, 
showing that this gain, as in the LOCC case, has a global maximum, and also decreases with larger systems sizes. 

In answering the question about existence of self-catalysis, we obtained many results, not only on deterministic 
and probabilistic self-catalysis, but also on ordinary catalysis. About the typicality of self-catalysis, our data support 
two conjectures: under LOCC, the probability of finding a self-catalytic reaction increases monotonically attaining a 
global maximum for a dimension about 20. We believe that the origin of this overall non-monotonic behaviour is due 
the competition between the sizes of each Schmidt vectors and the measure eoneentration phenomenon. On the other 
hand, the data suggests that under SLOCC the probability of finding a self-catalytic reaction increases monotonically 
with the dimension. 
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FIG. 7: Self-catalytic probability gain and direct conversion rates for randomly (Haar) chosen pairs of incomparable 
Schmidt vectors with size n = 45. The colour represents the number of pairs per pixel. 


In a sense, we estimated numerically the volumes of the sets of pairs of Schmidt vectors where the phenomena 
take place, but it was not possible to characterize completely the asymptotic behavior with the vectors sizes. In fact, 
some of our numerical results have a reasonable dependence on the way we sort the random Schmidt vectors. That is 
something to be explored elsewhere. Perhaps a better understanding of the geometry of the sets involved could help 
in that analysis. 
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